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Abstract 

We show how the resummation for time dependent quantities at high 
temperature can be performed with an effective classical theory. As 
an application we demonstrate that the leading term in the shear 
viscosity, which is related to the /9<j,2$2 spectral function can be calcu¬ 
lated classically, either using classical linear response theory or from 
the classical correlation function. The classical result depends ex¬ 
plicitly on the cutoff, and the choice A ~ T reproduces the known 
quantum result. 
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High temperature quantum field theories are known to suffer from a num¬ 
ber of IR problems which make the direct application of perturbation theory 
unreliable. Part of these problems can be cured by integrating out the hard 
thermal modes. For the static quantities this yields dimensional reduction 
[Q], for the nonstatic ones, depending on the way we do it, we arrive at the 
HTL effective action 0] or an effective classical field theory In par¬ 

ticular it was shown that in the <I>^ theory the self-energy can be calculated 
classically, and also the first quantum correction can be reproduced from the 
effective theory [|, 

The transport coefficients can be calculated from the microscopic the¬ 
ory using linear response theory |^. The shear viscosity gets the dominant 
contribution from the spectral function of the operator $^( 0 :) H 


T] = lim 

Po.P^O 


(po, p) 

Po 


= lim 

po,P^O 


Disc ([$ 2 (x), <1)2(0)]) 


Po 
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One can use a quantum theory with effective (resummed) spectral functions 
to compute this quantity 1,0,0]- The goal of this paper is to show that the 
same results can be obtained from the classical theory as well, which provides 
a simple calculational possibility as well as a feasible numerical framework. 

In the followings we shortly recall (cf. Ref.s i,|,0) and generalize, how 
one can use the classical theory to perform resummation for time dependent 
quantities. We will then apply the formalism to the viscosity. 

Resummation with Classical Theory 

Dimensional reduction, which yields an extremely powerful method to 
compute static quantities at high temperatures uses a very general, renor¬ 
malization group (RG) inspired technique to get rid of IR divergencies: it 
identifies the most IR sensitive degrees of freedom, separates them and inte¬ 
grates over the remaining ones. In this form it can be generalized to develop 
a resummation method also for the nonstatic quantities. To this end we have 
first to find the IR sensitive degrees of freedoms, with other words the source 
of the IR divergencies of the Feynman diagrams. The diagrams are generated 
by the generating functional; for the theory at finite temperature it reads 

0 
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where 


G{k,t) = — (cos cukt (3) 
2uJk ^ > ujk 

is the propagator, n{uj) = {e^‘^ — 1)“^ the Bose-Einstein distribution, c is a 
real time contour (eg. the Keldysh contour) and In the IR 

regime where |k| T the vacuum part behaves as ~ l/cn, the matter part 
is ~ T/u'^, because n{uj) = T/ui + 0{1). For massless fields the latter is 
quadratically, the former just linearly diverges in the IR. Let us denote by 
Gir the IR sensitive part of the propagator; it is defined asymptotically, its 
general form reads 


G'//?(k,t) = ^ (1 + C>(/3cc;)) coscukt. (4) 

Let us moreover denote the IR regnlarized propagator by G = G — Gir. 

The dimensional redaction technique suggests that we should rearrange 
the perturbation theory and postpone the calculation with the most singular 
propagator. We can try to represent the IR propagator by a Gaussian path 
integral, like in the RG, and then interchange the order of the operations. 
Symbolically, if 

J d(j) ~ 

then the effective theory is given by 
Z[j] = 

The perturbation theory generated by the interactions is now IR finite, and 
provides an action for its background field which will be the kernel of the 
subsequent path integral. In the standard RG we use the formula to lower 
the cutoff by considering G{k) = 0(A > \k\ > A — AA)G(A;), but in this way 
we can introdnce completely new type of degrees of freedom. 

In our case the very special form of the time dependence of the IR prop¬ 
agator allows a 3D representation. In fact we can write 
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where 


Ho[(j),Tr] = 
<1(0, Tt) = 


d^k 


(27r)3 2K^ 

/ r^k 

Ck 0k + - TTk 


TT-kTTk + Ci;k0-k0k 


(27r 


CJk 
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and, in the leading order 


Kk = T, Ck{t) = coscukt, ? 7 k(t) = sinwkt. (9) 


Since the IR propagator is dehned only in the leading order (cf. (|D), we can 
freely choose the snbleading parts of K, ( and r]. The different choices of the 
snbleading parts lead to different 3D theories, corresponding to the different 
prescriptions of the papers 1 S 0 - We will use here 


Rk = r0(A-(nk) (Ck(t),hk(t)) = 


(cos(a;kt),sin(cnkt)), if t G Ci ,2 
(1, 0), if f G C 3 , 


( 10 ) 

with the cutoff A ~ T which describes the validity range of the classical ap¬ 
proximation [0. The cutoff is useful also to suppress the spatial nonlocalities 
in the effective action [|^], and to ensure the quantum decoherence | 0 ]. 

To simplify notations we can introduce the effective generating functional 




u 
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where the normalizing factor Af assures Z[ 0 ] = 1 and thus cancels vacuum 
diagrams. It is background dependent, and we will include it into the Hamil¬ 
tonian 


Z|j] = / 


( 12 ) 


where H = Hq + In A/”. Because of the special form of the background held 
on the Matsubara contour (pH]) the effective Hamiltonian can be computed 


H 


r 


dim.red I 


+ 


2 


d^k 

(27r)3 


k^ki 


( 13 ) 


where T^jm.redM is the effective action of the dimensional reduction. 
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The background field, beeing free, can be easily generated by the initial 
conditions 

= — Jdt (14) 

where 

J(t) = (al'(to) + ^{to)dt)5c{t - to) (15) 

is the current localized at the initial time. This yields 

= (16) 


where J = j + J. The last term is current-independent, so it is canceled by 
the normalization. What remains is that Z[j,^\ = Z[J, 0], the background 
dependence can be absorbed into a redefined current. 

We can also perform a Legendre transformation with respect to this cnr- 
rent. The resulting effective action r[(y9] is background independent, it can be 
calculated using the ordinary real time perturbation theory with the propa¬ 
gator G. The IPI vertex functions can be obtained by differentiating T with 
respect to (p, and take it at the physical point ^Pphys which corresponds j = 0 


r(”)(a;i, ...Xn) 


<5fb] 


6ip{xi). ..6ip{Xn) 


—‘^phys 


(17) 


j = 0 means J = J. The Dirac deltas in J at the initial time set the initial 
conditions for the time evolntion. Therefore 




phys 


= 0 , 


^physi^o) (j) 
^t^physi^^o) 


(18) 


That is, while the averaging over the initial conditions corresponds to the 
local effective action, the time evolntion is governed by the time dependent 
effective action. The difference comes from the time-nonlocal loops 0, a 
gennine quantnm effect. 


Shear Viscosity in Scalar Field Theories 

As an application we can calculate the shear viscosity in theory. First 
we recall the quantum result R OT, then the different classical approaches. 
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The quantum result The calculation is not too involved even in the quan¬ 
tum case. The retarded Green functions have to be computed by the rule 


G%{x) = - G^2 b{x). (19) 

In the Feynman diagrams we have to use the propagators 

iG^\k) = iG^{k) + iG^^{k), iG^^{k) = n{ko)p{k), 

iG‘^‘^{k) = —iG^{k) + iG‘^^{k), iG'^^{k) = (1 -|- n{ko))p{k), 

( 20 ) 


where iG^(k, t) = 0(^)p(k, t)- We have performed here a self-consistent 
resummation and used interacting spectral function instead of the free one 
(cf. I, [im). This yields (including the symmetry factor 2) 


iG 


R 

(J)2<J)2 


(pG) = 2 


d^k 

(27r)3 


0(t) [p(k, t)p{p - k, t) + 2p(k, t)iG<{k, t)]. (21) 


Performing a Fourier transformation in time 

d^k du du' p(k,uj)p{p — k,uj') 


G^2^2{p) — 2 


(1 -1- n[uj) n{uj'y). (22) 


(27r)3 27r 27r pq — uj — uj' + is 
The discontinuity can be simply obtained using Disc(x -1- ie)~^ = —27ii6{x) 

r d'^k 

p$2$2(p) = 2 j p{k)p{p -k){l + n{ko) + n{po - ko)). (23) 

Using the identity 

1 -I- n{ix) + n{uj') = (1 — e“^*'‘^+‘^'^)(l -f n(a;))(l -|- n{uj')) (24) 

we get back the previous results i0- 


The classical approach As it is proven before, the same quantum result 
can be obtained using a two-step method, where in the hrst step we compute 
the effective operator only, with the IR stable propagator and in the back¬ 
ground of the IR helds. A typical diagram contributing to the classical (ie. 
without quantum loops) part for the retarded Greens function is shown on 
Fig. g. In formula 
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Figure 1: A typical classical contribution to the effective operator. 


Gf2<j,2^e//(^) = 4<F(a:) G^^{x, 0; $) <F(0), 


R 


(25) 


where 0; <F) is the retarded Greens function (cf. S0)- 

The same result can be obtained from the classical linear response theory. 
To see it let us add a current term to an action 


S(j) = S + 


(26) 


where / is an arbitrary function. We want to examine the linear response of 
the function g{^) to this change 


Gf(x-x') = 




6j{x') 


j=0 


= g'{<l>{0,x))Gf^{x - x'). 


We can get from the modihed equations of motion 


0 = 


SS{j) 


6S 


5<F(x) 5<F(x) 

by differentiating with respect to j: 
f 6S 




affix'-y) = -f my)) 6{x-y). 

) 


Its solution is simply 

G%(x - y) = Gi^{x - y) rmy)), 

SO hnally 

Gfg{x - x') = g'{^{x)) G^m - x') f{^{x')). 

For / — ^ we reproduce 


(27) 

(28) 

(29) 

(30) 

(31) 
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After the 3D integration the backgronnd lines are closed with 


{4(A^)4(«)> = (27T)‘S(k + 5) !G3d(*:). 


(32) 


The backgronnd helds, however, also follow a nontrivial, nonlinear time evo¬ 
lution, and the complete calculation cannot be performed in its generality. 
We will make an approximation similar to the one in the quantum theory: 
we join all the background lines on one propagator (i.e. sum up the the self¬ 
energy diagrams) and work further with these effective propagators. This 
yields the spectral representations 



(33) 


where we should use the complete spectral functions instead of the free ones. 
The retarded Greens function reads 



(34) 


which yields the discontinuity (note the symmetry fco ^ po ~ ^o) 



(35) 


There is another use of the approximation (|3^) , because it provides a 


simple way to extract the spectral function (cf. []^ for the quantum case) 


Pab(p) = PvoiGABfiDip)- 

In our case we have to compute the classical expectation value 


(36) 



(37) 


and, after Fourier transformation 



(38) 


which is indeed identical with (|3^. 
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Let us finally compute the viscosity from the classical theory. We approx¬ 
imate the spectral function by 


p{k) 


4fco7k 


UJ- 


2\2 


+ 16*:o7k’ 


(39) 


which yields in the narrow width approximation B0 

Ti5{kl - ujI) 


p{ky 


(^k7k 


(40) 


To be consistent we have to use the classical value of the damping rate. Its 
parametric form has been given by |B, it is equivalent to a high tempera¬ 
ture approximation of the parametric quantum result ||TT| , p!5[| . Therefore 

we can extract the on shell classical damping rate from the corresponding 
quantum expression [T^ performing a high temperature expansion 


7k = 


ISSbTTCUk 


1 ® (l 


(41) 


where L 2 = —/Q^ln(l — t)/tdt is the Spence function. After interchanging 
the order of integrations we can write 




ln(l -h t) 
t^/i. 1 



The integral vanishes for |k| =0; for large |k| we hnd 


7k(4.^kI |k|3>mT 


47oa;o [l + O i^ln — 
\ \ k ttit 



From (P, 


(^) and (|40|) we then obtain 


V 


p^o pq J (27r)4 fcg 

768 r 70^0 

A^Ttt J (k^ -|- mj.)T2 yi^cuk 



d^k 1 

(27r)3 7kW^ 


(42) 


(43) 


(44) 
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Since 7kCo’k is bounded, rj will have a logarithmic divergence, its coefficient 
can be extracted from the large momentum behaviour of 7ki-i^k (cf. (^31) ) 


rj 


192 

X^Tti 


, A 

In- + const. 

rriT 


(45) 


With A T we get back the leading term of the quantum result . Using 
the complete expression (El we hnd const= 2.1552. 

This divergence cannot be canceled by any local counterterm in the La- 
grangian. It is the consequence of having a composite operator which needs 
renormalization even in the classical case. The complete result, of course 
is independent on this auxiliary cutoff, the UV integration should carry the 
appropriate counterterms. 

Conclusion and Outlook 

We have summarized, how an effective classical theory can solve the re¬ 
summation problems in high temperature quantum held theories. For static 
quantities this classical theory is equivalent to dimensional reduction, the 
time evolution is governed by the effective quantum action. 

With this effective classical theory we could reproduce in the model 
the quantum result for the shear viscosity. It could be computed either by 
using the dehnition (from the retarded Greens function) or directly from the 
3D expectation value of The usual classical theory gives 

a cutoff dependent result even after a proper renormalization. This cutoff 
dependence has to vanish if we calculate also the UV contributions to the 
effective operator. 

The viscosity in this form is not complete, as shown in H- The correc¬ 
tion terms (ladder diagrams) can be of the same order as the leading one. 
Summing them up is far from beeing trivial. Since, however, the important 
1/A^ behaviour of the viscosity is essentially classical, one can try to perform 
the summation of the ladder diagrams in this approach, which may be easier 
than in the quantum theory. This is a task for future studies. 

The classical theory, on the other hand, can be simulated on computers. 
Since the viscosity is directly proportional to the <I>^ two-point function, 
it is a relatively simply accessible quantity. The simulations have to be 
performed with a hnite cutoff of the order of the temperature, or with the 
proper renormalization factor stemming from the UV integration. 
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